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Representation of ice clouds in radiative transfer simulations is subject to uncer- 
tainties associated with the shapes and sizes of ice crystals within cirrus clouds. 
In this study, we examined several ice cloud models consisting of smooth, rough- 
ened, homogeneous and inhomogeneous hexagonal ice crystals with various aspect 
ratios. The sensitivity of the bulk scattering properties and solar reflectances 
of cirrus clouds to specific ice cloud models is investigated using the improved 
geometric optics method (IGOM) and the discrete ordinates radiative transfer 
(DISORT) model. The ice crystal habit fractions in the ice cloud model may 
significantly affect the simulations of cloud reflectances. A new algorithm was 
developed to help determine an appropriate ice cloud model for application to the 
satellite-based retrieval of ice cloud properties. The ice cloud particle size retrieved 
from Moderate Resolution Imaging Spectroradiometer (MODIS) data, collocated 
with Multi-angle Imaging Spectroradiometer (MISR) observations, is used to infer 
the optical thicknesses of ice clouds for nine MISR viewing angles. The relative 
differences between view-dependent cloud optical thickness and the averaged value 
over the nine MISR viewing angles can vary from —0.5 to 0.5 and are used to evalu- 
ate the ice cloud models. In the case for 2 July 2009, the ice cloud model with mixed 
ice crystal habits is the best fit to the observations (the root mean square (RMS) 
error of cloud optical thickness reaches 0.365). This ice cloud model also produces 
consistent cloud property retrievals for the nine MISR viewing configurations 
within the measurement uncertainties. 


1. Introduction 

Field campaigns such as the First International Satellite Cloud Climatology Project 
(ISCCP) Regional Experiments (FIRE) (Starr 1987), Tropical Rainfall Measurement 
Mission (TRMM) (Silva Dias et al. 2002, Yuter et al. 2005) and Cirrus Regional Study 
of Tropical Anvils and Cirrus Layers - Florida Area Cirrus Experiment (CRYSTAL- 
FACE) (Jensen et al. 2004) have demonstrated that ice crystals within cirrus clouds 
are normally defective and irregular in shape (Heymsfield et al. 1990, Arnott et al. 
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1994, Korolev et al 1999, Heymsfield et al. 2002, Nasiri et al 2002, Heymsfield and 
Miloshevich 2003, Baumgardner et al 2005, Chepfer et al 2005). The great natural 
variability of cloud ice crystal habits causes large uncertainties in the computation 
of the transfer of electromagnetic radiation involving cirrus clouds and has drawn 
a great deal of attention, as evidenced by the large volume of research devoted to 
the topic (Macke 1993, Yang and Liou 1998, Masuda et al 2002, Mcfarquhar et al 
2002, Baran et al 2005, Lawless et al 2006, Wendisch et al 2007, Yang and Fu 
2009). To improve the treatment of the bulk radiative properties of ice clouds in vari- 
ous remote-sensing and climate modelling applications, the applied optics community 
needs to first fully understand and provide accurate scattering and absorption charac- 
teristics of a monodispersed ensemble of ice crystals, since this represents the simplest 
arrangement of ice crystals. 

A fundamental hindrance to simulating the light scattering of ice crystals is the dif- 
ficulty in obtaining realistic representations of natural ice crystals arising primarily 
from their complex particle shapes. Individual ice crystal structure is affected by a 
number of factors including the temperature, humidity and pressure at which ice crys- 
tals form (Kobayashi 1961, Kajikawa and Heymsfield 1989, Bailey and Hallett 2004, 
Baker and Lawson 2006, Lawson et al 2006). Several physical processes can further 
modify the appearance of ice crystals during their vertical motion within clouds. Ice 
crystals can become unidentifiable aggregates when they fall through the cloud layer 
and collide with supercooled water droplets or other ice crystals. Additionally, the 
surfaces of ice crystals are generally not perfectly smooth because of collisions, coa- 
lescence and other fragmenting processes. Although the exact mechanism for forming 
various ice crystal habits is not well understood, hexagonal columns (HCs) and plates 
(HPs) are known to be the important habits of ice crystals commonly found in cirrus 
clouds. This has been confirmed by both in situ measurements of cirrus clouds and 
the optical phenomena that 22° halos, caused by sunlight refracted by hexagonal ice 
crystals, are usually seen when the sky is overspread with thin cirrus clouds (Coleman 
and Liou 1981, Sassen et al 1994). 

Based on the ice crystal habit and size distributions obtained from in situ measure- 
ments and the work conducted to simulate light scattering by nonspherical particles 
(Cai and Liou 1982, Takano and Liou 1995, Macke and Mishchenko 1996, Yang and 
Liou 1996a, b, 1997, 1998, Baran et al 2001b), ice crystals with hexagonal structures 
are often assumed in the retrievals of cirrus cloud microphysical and optical proper- 
ties from satellite observations (Minnis et al 1993, King et al 1996, Minnis et al 1998, 
Chepfer et al 2001, King et al 2004, Minnis et al 2004, Li et al 2005, Minnis et al 
2011). A mixture of more complicated ice crystal shapes has also been suggested for 
ice cloud optical models (Baum et al 2000, 2005a, b, 2007, Xie et al 201 1). The quality 
of the satellite retrieval depends on the accuracies of the scattering properties used to 
represent natural hexagonal ice crystals in cirrus clouds and on the resulting simulated 
cirrus cloud solar reflectances. 

Efforts have been made to study the sensitivity of certain morphological parame- 
ters of hexagonal ice crystals to the single-scattering properties of ice crystals. Yang 
and Liou (1998) modelled the ice crystal surface roughness by using the slopes of 
the facets to define the degree of surface roughness. The radiation scattering by the 
rough-surfaced hexagonal ice crystals has been investigated using both rigorous and 
simplified algorithms based on the ray-tracing technique (Yang et al 2008b). The pri- 
mary effect of ice crystal surface roughness, relative to using smooth crystals, is to 
increase the solar reflectances in visible regions at most angles and to decrease the 
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retrieved ice cloud optical thicknesses (Yang et al. 2008a). Xie et al. (2009) investigated 
the scattering properties of inhomogeneous ice crystals based on observations, in 
which atmospheric ice crystals were collected near the surface and studied with a 
binocular microscope (Tape 1994). Inhomogeneous ice crystals, containing spherical 
or spheroidal air bubbles, were confirmed to have the effect of reducing the ice cloud 
optical thicknesses when they were used in the satellite-based retrieval of ice cloud 
properties. In addition, air bubbles within ice crystals may substantially reduce the ice 
crystal volumes and affect the retrieval of ice cloud particle sizes. 

Quantifying the morphological parameters (i.e. degrees of ice crystal surface rough- 
ness and inhomogeneity) and habit distributions of the ice crystals used in the 
operational retrieval of cloud properties would help minimize the effect of the spa- 
tial and temporal limitations of the in situ and surface observations of atmospheric ice 
crystals. In other words, a robust algorithm is needed to quickly and accurately find 
the most appropriate ice crystal configurations to minimize the uncertainties intro- 
duced by the use of any particular ice cloud model. Doutriaux-Boucher et al (2000) 
derived ice cloud optical thicknesses based on the Polarization and Directionality of 
the Earth’s Reflectances (POLDER) observations of radiances in 14 different view- 
ing directions. Cloud spherical albedos were computed using the retrieved optical 
thicknesses over ‘super-pixels’ consisting of 9 x 9 elementary pixels, each of which 
corresponds to a spatial resolution of approximately 60 km x 60 km. The error caused 
by the use of ice crystal habits and their scattering properties was inferred through the 
departure of the spherical albedos derived from their averaged values in conjunction 
with the 14 view angles. This algorithm was followed by studies aimed at improving the 
scattering properties of cirrus clouds (Labonnote et al 2001, Baran and Labonnote 
2006) or retrieving ice crystal habits (Baran et al 2001a, Chepfer et al 2001, 2002, 
Sun et al 2006) using POLDER multi-angle radiance measurements. McFarlane 
et al (2005) developed a method for retrieving cloud thermodynamic phase and ice 
crystal habit using Multi-angle Imaging Spectroradiometer (MISR) and Moderate 
Resolution Imaging Spectroradiometer (MODIS) data. They assessed the accuracy of 
ice crystal habits by studying the relative deviation of the computed cloud reflectances 
and those measured by nine MISR cameras and averaged over a 10 km x 10 km box at 
the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site. 
The retrieval of ice crystal habits was reasonably consistent with in situ observations 
of cirrus clouds. 

The differences in retrieved cloud optical properties caused by using measurements 
taken at different satellite viewing directions can arise from other factors besides 
ice cloud habits. Theoretical studies (Davies 1984, Loeb and Davies 1997, Loeb and 
Coakley 1998, Iwabuchi and Hayasaka 2002) have indicated that the horizontal inho- 
mogeneity of clouds or the 3D radiative effect may lead ID retrievals to under- or 
overestimate cloud reflectances and optical thicknesses due to cloud-side shadowing or 
illumination. However, because of its complexity, this effect has typically been ignored 
in current operational retrieval of cloud properties. Moreover, the biases in cloud prop- 
erty retrievals introduced from viewing geometries can be partially eliminated when 
the results are averaged over several pixels (Oreopoulos and Davies 1998, Kato and 
Marshak 2009). Thus, the improvement of pixel-by-pixel retrievals of cloud properties 
is more meaningful if the errors are investigated using similar spatial resolutions. 

In this study, we develop an ice cloud reflectance model using simple, but appro- 
priate, ice crystal habit and size distributions to minimize directional reflectance 
uncertainties employing matched data from the MODIS and MISR on board the 
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National Aeronautics and Space Administration’s (NASA’s) Terra spacecraft. Given 
all of the uncertainties of the operationally retrieved ice cloud properties, the resulting 
ice cloud model exhibits excellent consistency with the measurements from the nine 
MISR views. 

In §2, we introduce the study of ice crystal habits and their sensitivity to the 
scattering properties and solar reflectances of cirrus clouds. In §3, we describe the 
MODIS and MISR data used in the retrieval of ice cloud properties and develop 
a methodology to determine an optimal ice cloud model for routine retrieval of ice 
cloud properties using reflected solar radiation. In §4, we analyse the potential errors 
in retrieved ice cloud optical thicknesses in a case in which MODIS data are not 
employed in the retrieval. In §5, the errors caused by ice cloud models and their scatter- 
ing properties are addressed through analyses of viewing angle dependence of retrieved 
ice cloud optical thickness, and an ice cloud model with mixed ice crystal habits is 
selected from 200 models. Conclusions are provided in §6. 


2. Ice cloud models and their sensitivity to the simulations of scattering properties 
and solar reflectances by cirrus clouds 

2.1 Ice cloud models 

Nearly all atmospheric ice crystals develop on the fundamental hexagonal matrix that 
results from the near tetrahedral bonding angles of the water molecule. Hexagonal ice 
crystals, in one form or another, are one of the most important ice crystal shapes 
within cirrus clouds. In this study, only HCs and HPs are assumed in developing 
an algorithm to determine optimal ice cloud models for satellite-based remote sens- 
ing of cirrus clouds. Other ice crystal habits such as bullet rosettes, droxtals, hollow 
columns and aggregates, often observed in clouds, could also be used for the ice cloud 
models to, perhaps, further improve their performance in retrieving ice cloud proper- 
ties (Baum et al. 2003, 2005a, b). The habit fractions of ice crystals may have regional 
preferences and are strongly related to cloud temperature. This aspect may deserve 
future study but is beyond the aim of this work. 

Hexagonal ice crystal geometries are normally specified by their aspect ratios, 
a = 2 a/ L, where a and L are the semi-width and length of the ice particles, respec- 
tively. The aspect ratios of HCs used in this study are consistent with those of Mitchell 
and Arnott (1994) and can be formulated as follows: 


0.7 for D < 100 (xm 
6.96L~ 05 for D > 100 |xm ’ 


( 1 ) 


where D represents the dimensions of the ice crystals and equals L for HCs. For HPs, 
the aspect ratios are determined based on in situ measurements of cirrus clouds (Auer 
and Veal 1970, Pruppacher and Klett 1980) and can be given by 


a = 


1 for D < 4 |xm 

0.2227a + 1.5547 f ° r 
0.8038«°- 526 for D>10|xm 


( 2 ) 


where D = 2a. 

Ice crystals with roughened surfaces or internal air bubbles are also considered here 
because they retain a simple hexagonal shape while significantly altering the optical 
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properties. As in Yang et al. (2008b), the root mean square (RMS) tilt, cr, of the 
roughened facets on the ice crystal surfaces is used to specify the degree of roughness. 
Thus, the surface condition varies from smooth to severely rough as a increases from 
0 to 1 . The geometries of inhomogeneous ice crystals have been described in detail by 
Xie et al. (2009) and will not be restated here in the interest of brevity. The internal 
air bubbles of the inhomogeneous ice crystals are simulated using spheroidal shapes 
that satisfy 



(3) 


where x, y and z represent the axes of x, y and z, respectively, and n , U and rj are the 
air bubble semi-axes along the x, y and z directions, respectively. Here, n, ty and r-? 
depend on the widths and lengths of the ice crystals as follows: 


D = Cbubbie^? 

(4) 

— Cbubbie^? 

(5) 

r 3 = Cbubbie D, 

(6) 


where Cbubbie is a constant. 

For the ice cloud models, the ice crystal surface roughness condition and internal 
air bubbles are used for both the HC and HP cases. The geometric parameters a and 
Cbubbie are set to be 1 .0 and 0.5, respectively, indicating a case with severely roughened 
ice crystals and a case with medium-sized air bubbles. Thus, the scattering proper- 
ties of cirrus clouds are represented by ice cloud models consisting of six types of 
ice crystals including HC, surface-roughened column (RC), inhomogeneous column 
(IC), HP, surface-roughened plate (RP) and inhomogeneous plate (IP). Figure 1 com- 
pares the projected area and volume of these ice crystals as functions of ice crystal 
dimension. The degree of ice crystal surface roughness has not been considered in the 
computation of the projected areas and volumes of the ice crystals. As can be seen 
in figures 1(a) and (b), the projected areas of HP are much larger than HC when the 
dimensions are greater than 250 qm. However, the volume of HC is greater than HP 
in the region of 0-500 |xm, as shown in figures 1(c) and (d). Note that the internal 
air bubbles within the ice crystals substantially reduce their volumes while having no 
effect on the projected areas. 

To quantitatively characterize the cirrus cloud particle sizes, we use the effective 
particle size from 


3 f V total (D)n(D)dD 

r\ J C/min / H\ 

I A total (D)n(D)dD 

^ T^min 

where n(D) is the ice particle size distribution, V ioia \{D) and A i0Vd \{D) are the total 
volumes and projected areas of the ice crystals for size D, respectively, and D max 
and D m in are the maximum and minimum values of D (Foot 1988). Following Baum 
et al. (2005a), the size distribution of ice cloud particles is represented by the gamma 
distribution given by 


Downloaded by [Texas A&M University Libraries and your student fees] at 08:16 22 December 201 1 


4224 


Y. Xie et al. 



Figure 1. Ice crystal projected area versus dimension: (a, b) for HC (black) and HP (green) 
models and (c, d) variation of ice crystal volume versus dimension for HC (black lines), 
IC (blue), HP (green) and IP (red). 


n(D) = No I) 1 ' exp 


b + /i + 0.67 

-^median 


D 


( 8 ) 


where No is the intercept, D ni edian is the median of the distribution of D, and /x and 
b are constants assumed to be 2.0 and 2.2, respectively. V loRi \{D) and A loVd \(D) in 
equation (7) are given as 

Ktota m = Chc Vh c(D) + Crc V rc (D) + Cic Vic(D) + 

Chp V hp (D) + C RP Frp(D) + C IP V lP (D), W 


^total(^) = Chc^Hc(^) + Crc^RcC^) + CicAic(D) + , , 

Chp^hp(^) + CrpHrp(D) + CipHip(D), 

where V HC (D), V RC (D), V ]C (D), V HP (D), V RP (D) and V lP (D) and A nc (D), A rc (D), 
A\c(D), ^hp!^), ^rp(D) and A\ P (D) are the volumes and projected areas, respectively, 
of the aforementioned ice crystal habits, denoted here by the subscripts; 6 'hc (^), 
Crc(D), Cic(D), Chp(D), Crp(D) and Cip(D) are the corresponding habit fractions 
of the ice crystals. Thus, the goal of this study was to determine the appropriate 
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coefficients Chc( D), Crc(D), C\q(D), Chp(D), Crp(Z>) or C\p(D) by using satellite 
observations on a global scale. 

2.2 Sensitivity of ice cloud models to the simulations of their scattering properties 

For single ice crystals, the light scattering properties are simulated by an improved 
geometric optics method (IGOM) (Yang and Liou 1996b). After simulation, the scat- 
tering properties including scattering phase functions, P\\(0), extinction efficiencies, 
Q e , and single-scattering albedos, co, of an ensemble of ice crystals can be derived on 
the basis of ice crystal habit and size distributions as follows: 


p Anax 
^ ^min 

Y J C i PnA&,D)QUD)A i (D) 

_ / 

n(D)dD 

J 

‘ Anax 
n 

^min 

E CiQUD)Ai(D) 

_ i 

n(D)dD 


r Anax 
A min 

E CiQ.mAfD) 

_ i 

n(D)dD 

f 

JD 

max 

nin 

1 1 
-M 

Ci 

b 

1 1 

n(D)dD 


r Anax 
'A D m [ n 

E QQ s AD)Ai(D) 

_ i 

n(D)dD 

pAnax 
** ^min 

E CiQ e j(D)Ai(D) 

- i 

n(D)dD 


( 11 ) 


( 12 ) 


(13) 


where 0 is scattering angle, P\\j{0) is the scattering phase function; Q e i and Q si are 
the extinction and scattering efficiencies, respectively, for each ice crystal habit, and 
the subscripts i represent HC, RC, IC, HP, RP or IP. 

Figure 2 compares the scattering phase functions of the ice cloud models for 
D e = 50 pm. To investigate the sensitivity of the phase function to ice cloud modes, 
each ice cloud model consists of only one ice crystal habit. Thus five of the habit 
fractions are 0 for each ice cloud model. From figure 2, HPs are seen to generally 
have larger forward scattering than HCs, although both ice cloud models have the 
same cloud particle size. Most HPs have larger projected areas than HCs for the same 
dimensions (see figure 1). Therefore, the diffraction of HPs is more significant than the 
scattering of light by the ice crystals. At a wavelength of k = 0.866 pm (see figure 2(a)), 
the phase functions of HCs and HPs have strong forward and backward scattering 
and distinct halo peaks at the scattering angles near 22° and 46°. When surface rough- 
ness is introduced (i.e. RC and RP), the phase function halo features, caused by solar 
beams refracted by hexagonal ice crystals, completely vanish. Moreover, the forward 
and backward scattering of HCs and HPs are reduced due to light beam spreading 
in the refraction processes of the ice crystals (Yang et al. 2008b). As discussed by Xie 
et al. (2009), the inhomogeneity effect on the phase functions can also be seen for 
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Figure 2. Scattering phase functions for HC (black), RC (blue), IC (cyan), HP (green), RP 
(yellow) and IP (red) models composed of randomly oriented ice crystals. 


single ICs and IPs. However, the volumes of ICs and IPs are smaller than HCs and 
HPs because the volumes of air bubbles are subtracted. To achieve the same cloud 
effective particle sizes, inhomogeneous ice crystals in the ice cloud models have larger 
dimensions and larger aspect ratios than the homogeneous ones. Thus, the forward 
scattering of the IC and IP models is stronger than that for the HC and HP models, 
as shown in figure 2(b). Figures 2(c) and ( d ) compare the six phase functions of the 
ice cloud models at the wavelength of A. = 2.13 |xm. At the near-infrared (NIR) wave- 
length, the smoothing of the scattering phase functions is shown for ice cloud models 
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consisting of RCs, RPs, ICs and IPs. At the NIR wavelength, the smoothing of the 
scattering phase functions for the RC, RP, IC and IP models relative to the homoge- 
neous models is similar to that seen at 0.866 pm. As X increases from 0.866 to 2.13 
|xm, the phase functions of the ice cloud models decrease in the forward scattering 
directions, because the ice crystal size diminishes relative to the wavelength. 

2.3 Sensitivity of solar reflectance to cirrus ice particle habit 

To retrieve ice cloud optical and microphysical properties, simulated cloud bidirec- 
tional reflectances are typically matched with satellite observations at solar wave- 
lengths. Because the simulated cloud reflectances are sensitive to the particle habits 
used in the models, the retrieved properties will vary with the ice crystal scattering 
model used. This sensitivity needs to be quantified to understand its impact on the 
retrievals. The satellite-observed cloud bidirectional reflectance, R, is defined by 


where F denotes the direct solar flux density at the top of the atmosphere (TOA), /u, 
and p are the cosine values of the solar and satellite viewing zenith angles, 9 o and 9, 
respectively, 4> is the relative azimuth angle, and /( p . </>) is the radiance observed by 
the satellite instruments. The relative azimuth angle 4> can be specified as follows: 


where cpo and <p are the solar and satellite viewing azimuth angles, respectively. 

Figure 3 shows the simulated cloud bidirectional reflectances using the six ice cloud 
models. The cloud reflectances were computed using discrete ordinates radiative trans- 
fer (DISORT) (Stamnes et al. 1988) assuming a homogenous ice cloud layer over a 
black surface, a cloud optical thickness of 1 and 9 o = 30°. In each panel of figure 3, 
9 increases from 0° to 90° along the radial direction, and <fi varies along the angular 
direction from 0° to 180°. At a wavelength of X = 0.866 pm. figure 3(a) clearly shows 
that the cloud reflectance is sensitive to viewing geometries. At 9 o = 30°, the maxi- 
mum values of the cloud reflectances are associated with the forward directions where 
9 is between 80° and 90°. In the backscattering directions, the cloud reflectance max- 
ima occur at 9 = 30° and are correlated to the phase function maxima at 0 = 180°. 
Clearly, the cloud reflectances in figure 3(a) are sensitive to the use of ice cloud models. 
For 9 = 0°-45°, a range frequently used for satellite-based retrievals, the HP, RP and 
IP reflectances are much smaller than their columnar counterparts. This is probably 
due to the stronger forward scattering by plates (figure 2). The relative maxima at 
0 = 180° in the HC and HP phase functions (figure 2) yield the nearly pinpoint 
backscattering reflectance maxima seen in both figure 3(«) and (b). However, these 
maxima are significantly reduced for the inhomogeneous ice crystals and absent for 
the rough crystal reflectances as a result of their smoother phase functions. The sensi- 
tivity of cloud reflectances to the choice of ice cloud model is evident in the simulations 
at both wavelengths. Therefore, the ice crystal habit fractions in the ice cloud models 
will significantly affect the simulations of cloud reflectances and the retrieval of ice 
cloud properties. 



(14) 


| cp — <pol for 0° < \cp — <po\ < 180° 

360° — | cp — cpo\ for | cp — cpo\ > 180° ’ 


(15) 
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000 °- 10 °- 20 >0.30 o.OO 0.05 0.10 0.15 >0.20 

Figure 3. Cloud bidirectional reflectances computed for the HC, RC, IC, HP, RP and IP ice 
cloud models for 6 0 = 30°, D e = 50 pm and (a) X — 0.866 pm and (b) X — lA'i pm. 


3. Data and methodology 

3.1 Terra MOD IS and MISR data 

Terra, launched on 18 December 1999, is a scientific research satellite component of 
the Earth Observing System (EOS). The MODIS instrument is designed to monitor 
the state of the Earth’s environment including the radiation processes occurring in the 
atmosphere and at the Earth’s surface. MODIS has 36 bands ranging in wavelength 
from 0.414 to 14.235 pm (Ardanuy et al. 1991). In this study, the MODIS bands 
centred at 0.86 and 2.13 pm are used to derive the ice cloud optical thickness and 
effective particle sizes over ocean because of the relatively small oceanic reflectance 
and atmospheric absorption at those wavelengths (Pinker and Laszlo 1992, Gatebe 
et al 2005). In the retrieval process, the viewing and illumination geometries and 
reflectances at the TOA are provided by the MODIS Level IB Calibrated Geolocation 
Data Set (MOD02) at a spatial resolution of 1 km x 1 km. The cloud thermodynamic 
phase is taken from the 1 km resolution cloud product (MOD06). The MOD02 and 
MOD06 data are stored in MODIS granules, each consisting of a 5-min data 
swath. 

The Terra MISR instrument is designed to simultaneously provide measurements of 
daytime shortwave radiances from nine distinct cameras that view the nadir direction 
and the forward and aftward local vertical directions (Diner et al 1989, 1998, 2002). 
MISR’s cameras are designated as Df, Cf, Bf, Af, An, Aa, Ba, Ca and Da, where n 
denotes the nadir direction and f and a indicate whether the camera points forward or 
aftward along the Terra ground track. Except for the nadir camera, D, C, B and A rep- 
resent the satellite viewing zenith angles of 70.5°, 60.0°, 45.6° and 26.1°, respectively. 
Each of the nine MISR cameras views the Earth in four spectral bands with central 
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wavelengths at 0.446 (blue band), 0.558 (green band), 0.672 (red band) and 0.866 pm 
(NIR band) (Diner et al. 2002). 

The Sun-synchronous Terra orbits have an exact 16-day repeating cycle. Because 
of this periodicity, MISR Level 1-2 data and data processing are divided into 
233 discrete paths for the entire mission period. Each MISR path is further divided 
into 180 blocks with block 1 starting in the Arctic and block 180 ending in the 
Antarctic. Each block has a fixed sub-orbital area and corresponds to a fixed lati- 
tude, but the block number is a non-uniform function of latitude. The MISR Level 
1B2 data (MI1B2T) and the retrieved reflecting level reference altitude (RLRA) 
values are used to develop the MISR Level 2 cloud product (MIL2TCAL), where 
the multi-angle reflectances from the nine cameras are registered at the height of 
the reflecting level (Diner et al. 1997, 1999). The 2.2 km resolution MIL2TCAL 
data are employed to derive ice cloud optical thicknesses using multi-angle cloud 
reflectance measurements. 

Figure 4 shows an RGB (red-green-blue) image of MODIS radiances over the cen- 
tral Pacific Ocean acquired at 2050 Coordinated Universal Time (UTC) on 2 July 
2009. The image was created using a MODIS Level IB granule, which corresponds 
to 1354 x 2030 1 km pixels. The white boxes are the simultaneous MISR observation 
domains for blocks 78-95 in Path 60. In figure 4, the MISR region can be seen over- 
lapped by the MODIS granule. By matching the MODIS data at 2.13 pm with the 
overlapping MISR data, it is possible to provide additional information, for example, 
particle size, to the MISR retrievals of ice cloud optical properties. 



Figure 4. MODIS RGB (R for band 1; G for band 4; B for band 3) image over the central 
Pacific Ocean, 2050 UTC, 2 July 2009. The white boxes indicate blocks 78-95 of MISR Path 60. 
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3.2 Algorithm to retrieve ice cloud optical thickness using MODIS and MISR data 

The first step in determining the optimal ice cloud model is the retrieval of the ice 
cloud optical and microphysical properties. The cloud reflectances are simulated in 
MODIS bands 2 and 7 and the MISR NIR band for the cloud optical thickness 
r = 0.05, 0.1, 0.2, 0.5, 1.0, 2.0, 4.0, 6.0, 8.0, 10.0, 12.0, 14.0, 16.0, 18.0, 20.0, 25.0, 30.0, 
35.0, 40.0 and 50.0. The cloud optical thicknesses at X = 0.866 |im (hereafter referred 
to as ro.866) are used as the reference optical thicknesses. The optical thicknesses at 
other wavelengths are related to ro.866 by 

r = ^ ro.866, (16) 

(ce,0.866 

where Q t represents the extinction efficiency of a given ice cloud model and 0 e ,o.866 
is Q q at A = 0.866 pm. The solar and satellite viewing zenith angles are constructed 
using /xo//x = 1.0, 0.95, . . ., 0.0, and the relative azimuth angles 0 are 0°, 2.5°, 5.0°, 
10.0°, . . ., 175.0°, 177.5° and 180.0°. 

Figure 5 shows a flow chart of the algorithm for retrieving ice cloud effective par- 
ticle sizes using MODIS data. The simulated cloud reflectances, Ro.866 and R 2 . 13 , 
are first interpolated in cloud optical thickness, particle size and solar and satellite 
viewing angles. For a cloud within a MODIS pixel, the thermodynamic phase is pro- 
vided by the MOD06. The effective particle size of an ice cloud is initially assumed 
to be 100 |xm. The value of ro.866 is retrieved by matching the MOD02 band 2 
reflectance to the reflectance simulated for each ice cloud model associated with D e = 
100 pm. A new effective particle size, D e 2 , of the ice cloud is then computed using 
the MOD02 band 7 reflectance, and 72.13 is determined with equation (16) using the 
retrieved ro.866- The final retrieval of the cloud particle size is equal to its initial guess 
(Dei) when |D e i - D e 2 | is negligible. Otherwise, the retrieval process is repeated, using 
an adjusted D e i, until the cloud effective particle size is retrieved or a specific stopping 
criterion is met. 

Figure 6 shows the cloud thermodynamic phases and D e retrieved using the HC 
model reflectances for the image in figure 4. The MODIS-retrieved cloud particle sizes 
within the MISR blocks can be used to retrieve the ice cloud optical thicknesses from 
MISR data. Figure 7 outlines the process for retrieving ro.866 at different satellite view- 
ing geometries. The retrieval is limited to MISR pixels over ocean overcast by ice 
clouds. The MODIS-retrieved cloud effective particle sizes, averaged over four con- 
tiguous MODIS pixels, are collocated with MISR pixels in 2.2 km x 2.2 km regions, 
when more than 70% of the MISR pixels are covered by the MODIS pixels. The 
4-pixel D e average is also compared with those for the surrounding 12 MODIS pixels 
to test the scene for uniformity. The cloud optical thickness retrieval is not performed 
when the D e difference to the surrounding pixels exceeds 20 pm. The bias in cloud 
optical thickness caused by using an inappropriate D e is discussed in §4. The collected 
D e is used in cooperation with the MISR NIR band reflectance of MIL2TCAL 
and the simulated Ro.866 to derive ro.866 for each MISR viewing geometry. 

4. MODIS retrieval dependence of cloud optical thicknesses retrieved by MISR 

The MODIS-retrieved cloud optical thicknesses are subject to uncertainties indepen- 
dent of the MODIS retrieval. As discussed in §3, cloud optical thicknesses and effective 
particle sizes are retrieved using cloud reflectances at visible and NIR wavelengths. 
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Figure 5. Flow chart of the algorithm for retrieving ice cloud effective particle sizes using 
MODIS data. 


In the retrieval of r, the derived cloud particle sizes are employed to match the model 
and observed reflectances. On the other hand, the retrieved cloud optical thicknesses 
can be used to compute the MODIS 2.13 pm reflectance to enhance the retrieval 
of D e . This process can be applied iteratively to improve both the r and If retrievals 
until the simulated cloud reflectances match the satellite observations. Therefore, 
biases in r may be present when only the MISR NIR band reflectances are used 
in the retrievals. 

To illustrate such errors, figure 8 (a) plots the lookup table reflectances for the 
MODIS band 7 (black lines) and the MISR NIR band (red lines) for the HC model 
for two relative azimuth angles 85° (figure 8(a)) and 177.5° (figure 8(b)) where /x o 
and tx are both set to be 0.85. The red curves (MISR NIR reflectances) were com- 
puted using D e = 50 pm. In figure 8(a), the retrieved cloud optical thicknesses are 
slightly overestimated when D s < 50 pm. For clouds consisting of larger ice crys- 
tals, the retrieved cloud optical thicknesses with MISR NIR band reflectances are 
significantly underestimated. The bias caused by the absence of the MODIS band 7 
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Figure 6. Cloud properties determined from radiances corresponding to the image in figure 4. 
(a) Thermodynamic phase and (b) ice cloud effective particle sizes derived using the HC ice 
cloud model. 


reflectances varies with the ice cloud models and solar and satellite viewing geome- 
tries of the lookup tables. The lookup table reflectances plotted in figure 8 (b) reveal 
that the errors in the cloud optical thicknesses using MISR data with D e = 50 (xm are 
more extreme when retrieved near the backscattering direction. 
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Figure 7. Flow chart of the algorithm for retrieving ice cloud optical thickness using matched 
MODIS and MISR data. 

Figure 9 shows the relative errors in retrieved cloud optical thickness (RE r ) using 
the six ice cloud models with D e = 50 |xm. These errors are defined as 

RE. = 1 V ™ 1SR -' - Ti 100%, (17) 

N > fir z < 

where r, is the cloud optical thickness used to simulate the MODIS band 7 and 
MISR NIR band reflectances, tmisr,/ is the optical thickness retrieved using the 
lookup table of MISR NIR band reflectances, with T representing the number of 
the collocated pixels, and N r is the total number of r, for each D e . Figures 9 (a)-(d) 
show RE r averaged over the scattering angle, 0, ranges 0°-45°, 45°-90°, 90°-135° 
and 135°-180°, respectively. In the region of forward scattering (see figure 9(a)), 
RE t is most sensitive to ice crystal habit and size. Fortunately, reflectances at 0 = 
0°-45° are infrequently used in satellite-based remote sensing. The RE r at 0 = 45°- 
90° is within ±10% for most of the models and particle sizes. For 0 = 9Q°-135°, 
the RE r values for the column models are nearly 0% for D t < 60 |xm and reach 
—20% for the largest crystals. The RE T values for the plate models in figure 9(c) 
are positive for D e = 0-50 and 120-180 |xm, indicating an overestimate of r 
for such clouds. For D e = 50-120 qm, RE T is slightly smaller than 0. In the 
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Reflectance (MISR NIR band) 

Figure 8. MODIS band 7 and MISR NIR band reflectances (black) for variable particle size 
and optical thickness and MISR NIR band reflectances (red) for D e — 50 |xm. All calculations 
performed for the HC model at /r 0 = /x = 0.85 and <f) — (a) 85°, ( b ) 177.5°. 

region of backscattering (figure 9(d)), RE r for the column models is negligible. 
However, the optical thicknesses of plates are significantly underestimated when 
D e = 50 |xm is used in the retrieval. 

From the results shown in figures 8 and 9, the average RE T values for the column 
models are within ±10% for the scattering angles related to satellite-based remote 
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Figure 9. Relative errors in r using HC, RC, IC, HR RP and IP models with D e = 50 |im and 
e = (a) 0°M5°, ( b ) 45°-90°, (c) 90°-135° and (d) 135°-180°. 


sensing. For plates, using an estimated cloud particle size may lead to substantial 
errors, especially for ice clouds with extremely small or large sizes. Thus, employing 
the MODIS retrieval of ice cloud effective particle sizes can efficiently remove the 
MISR retrieval biases that would otherwise result from the assumption of a constant 
particle size in the retrievals. 


5. Results 

Based on the algorithm described in §3, r is retrieved for the nine MISR viewing 
geometries for the matched MODIS and MISR pixels. The difference between the 
cloud optical thicknesses from the nine viewing angles and their average can be used 
to evaluate the accuracy of the retrieval. With the presence of other uncertainties in 
the retrieval of ice cloud optical properties, the most appropriate ice cloud model will 
give the smallest optical thickness difference (OTD) values and RMS error s of optical 
thickness and are formulated as follows: 
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OTD, v 


T iJ ~ <T i>j 


<Ti>j 


(18) 


s = 


\J 


Vp 


9Np 


EE 0Tr C 


7 = 1 i = 1 


(19) 


where i is the number of the MISR viewing geometries, j is the number of collocated 
pixels, t ij is the retrieved ice cloud optical thickness for each available MISR pixel, 
<t j>j is the averaged value of r ; / - over the nine MISR viewing geometries and Np is 
the number of available MISR pixels. 


5.1 Errors in retrieved optical thickness for matched MODIS and MISR granule 
pixels 

Table 1 lists the values of e for the matched MODIS and MISR data from the granules 
shown in figure 4 and those over the central Pacific Ocean on 19 February 2010 at 
2140 UTC. For the retrievals using data from the granule in figure 4, the smallest error 
is clearly that for the RC model. The IP model yields the best performance in retrieving 
r for the 1 9 February case. However, the retrieval errors for the HP, RP and RC models 
are not much greater than the IP minimum. Because of the natural variability of the 
cloud particle habits, the differences between the 2 July and 19 February results are 
not surprising. 

Figure 10 shows the histograms of D e and r retrieved using the HC model for 
the matched MODIS and MISR data for the granules taken on 2 July 2009 and 
19 February 2010. The distributions of D e and r for the data in figure 4 are centred 
at around 70 and 15 |xm, respectively. For the retrievals using data from 19 February 
2010, the particle sizes of the clouds are slightly smaller than 70 |im, and the maximum 
optical depth occurs around r = 3. 

Because the RC model error is small for both days and has the least error for both 
days combined, it would be the most appropriate model overall for these two cases. 
The IP model has the lowest error for the 19 February case: its overall error is ~19%, 
larger than that for the RC model. Because the IP model is the most appropriate for 
the 19 February case and the retrieved optical depths are, on average, considerably 
less than those retrieved for the 2 July case, it would appear that the IP or any of the 


Table 1 . RMS error, s, of ice cloud optical thickness retrieved from matched 
MODIS and MISR data for one 5-min MODIS granule using six different 

models. 


Ice cloud model 

2050 UTC, 2 July 2009 

2140 UTC, 2 July 2009 

HC 

0.2450 

0.198 

RC 

0.1760 

0.168 

IC 

0.2560 

0.195 

HP 

0.2576 

0.163 

RP 

0.2376 

0.166 

IP 

0.2446 

0.155 
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Figure 10. Flistograms of the ice cloud (a) effective particle size and ( b ) optical thickness 
retrieved using the HC model for the granules taken on 2 July 2009 and 19 February 2010. 


plate models would be best for optically thin clouds. However, the RC model error is 
only 8% greater than £(IP) for that case. Thus, the use of plates to optimally retrieve 
r for optically thin clouds is somewhat ambiguous. Further classification of the errors 
according to optical depth range would be required to determine whether the model 
errors vary with r . 

5.2 Examination of ice cloud models using 1-day matched MODIS and MISR data 

To further examine the potential of finding an ice cloud model that produces the 
least biased retrievals of ice cloud properties, models containing various mixtures 
of the six ice crystal models were tested using MODIS and MISR data on a global 
scale. Following the results above using the two granules of matched MODIS and 
MISR data, RC data are assumed to be the major component for the habit mix- 
tures. The habit fractions of the ice crystals in the ice cloud models are specified as 
follows: 


C RC = 0.4 + 0.6&, 


( 20 ) 
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C HC = 


(1 — Crc )^2 

Hi + £3 + §4 + Hi + He 


( 21 ) 


Ci C = 


(1 — Crc)§3 


Hi + $3 + H4 + £5 + He 


( 22 ) 


Chp = 


(1 — Crc)^4 


Hi + Hi + H4 + £5 + §6 


(23) 


Crp = 


( 1 — Crc)^5 


Hi + £3 + H\ + ?5 + He 


(24) 


Cip 


(1 — Crc)^6 

Hi + £3 + §4 + §5 + He 


(25) 


where £ 1 , £ 2 , Hi, H 4 , H 5 and £6 are independent random numbers uniformly distributed 
in [0, 1], 

Two hundred sets of habit fractions were initially examined using the two cases of 
2 July and 19 February. Table 2 lists the habit mixture models having the minimum 
RMS errors in r for Crc varying within the fractional intervals (0.4, 0.5), (0.5, 0.6), 
(0.6, 0.7), (0.7, 0.8), (0.8, 0.9) and (0.9, 1.0). The RMS errors of cloud optical thickness 
generally decrease with increasing Crc, which is associated with the overall low error 
(~0.172) of the RC model for the two cases combined. Including some fraction of 
other habits in the mixture only reduces the minimum two-case error by 0.6% (models 
Ml and M2) compared to the pure RC model. 

Among the ice cloud models in this study, those consisting of a single ice crystal 
habit and Ml, M3, M5, M7, M9 and Ml 1 are further tested using all daytime matched 
MODIS and MISR data taken over ocean on 2 July 2009. A total of 129 MODIS 
granules and 15 MISR paths are involved. Figure 1 1 shows the histograms of D e and 
r retrieved using the HC model applied to the global 2 July data set. Compared to 
figure 10(a), the distribution of the cloud particle sizes in figure 11(a) is populated 
with ice particles covering a broad range with a peak at 60 pm. For the retrieved 
cloud optical thickness, the peak of the distribution occurs near 4, with ~50% of the 
population having r < 15. The retrieved ice cloud properties in figure 1 1 are based on 
244 61 1 2.2 km MISR pixels, which indicates that more than 100 and 2000 retrievals 
contribute to any interval displayed in figure 11(a) and ( b ), respectively. The large 
number of samples from all over the globe confirms that the results for this 1-day 
data set are much more representative of ice clouds over the globe than the two cases 
of 2 July and 19 February. 

Figure 12 compares the ice cloud optical thicknesses derived using MISR measure- 
ments at nine different viewing angles. The retrievals of r for figure 12 were obtained 
using the RC model. The optical thickness for each MISR viewing angle i, <r>,, is 
specified by 
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Table 2. Ice cloud models and RMS errors of cloud optical thickness for the MODIS and MISR 
granules on 2 July 2009 and 19 February 2010. 


Ice cloud 
model 

Cue 

Crc 

Cic 

Chp 

Crp 

Cip 

£ 

Ml 

0.0001 

0.9974 

0.9 < C RC 
0.0014 

< 1.0 
0.0009 

0.0002 

0.0000 

0.171 

M2 

0.0023 

0.9828 

0.0034 

0.0067 

0.0003 

0.0045 

0.171 

M3 

0.0423 

0.8913 

0.8 < C RC 
0.0004 

< 0.9 
0.0115 

0.0063 

0.0482 

0.171 

M4 

0.0004 

0.8848 

0.0115 

0.0063 

0.0485 

0.0485 

0.172 

M5 

0.0353 

0.7963 

0.7 < C RC 
0.0390 

<0.8 

0.0598 

0.0158 

0.0538 

0.172 

M6 

0.0270 

0.7875 

0.0584 

0.0723 

0.0045 

0.0503 

0.172 

M7 

0.0589 

0.6987 

0.6 < C RC 
0.0361 

<0.7 

0.0710 

0.0374 

0.0979 

0.173 

M8 

0.0068 

0.6817 

0.0759 

0.0807 

0.0228 

0.1321 

0.173 

M9 

0.0701 

0.5217 

0.5 < C RC 
0.0191 

<0.6 

0.0779 

0.1609 

0.1503 

0.174 

M10 

0.0190 

0.5389 

0.0775 

0.1602 

0.1495 

0.0549 

0.174 

Mil 

0.0091 

0.4818 

0.4 < C RC 
0.0066 

<0.5 

0.0254 

0.2528 

0.2243 

0.174 

M12 

0.0057 

0.4026 

0.0221 

0.2196 

0.1948 

0.1552 

0.174 


<r)i = 


l 

Vpi 


l 

Np 2 


1 

Np3 


1 

JVp4 


1 

JVp5 


1 

Npe 


1 

Npi 


1 

Np 8 


Npi 

E Ti J 

7=1 

Npi 

E 

7=1 

Npj 

E T u 

7 = 1 

Npn 

E T ‘d 

7 = 1 

JVp5 

E r o- 

7 = 1 

Np& 

E r o- 

7 = 1 

Npi 

E r D- 

7 = 1 

JVp8 

E T D' 

7 = 1 


for (r) < 1 
for 1< (r) <2 
for 2< (r) < 3 
for 3< (r) < 4 
for 4 < (r) < 7 
for 7 < (r) < 15 
for 15 < (r) < 30 
for 30 < (r) 


(26) 
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Figure 11. Same as figure 10, except for 1-day data taken over ocean, 2 July 2009. 


where <r> is the averaged value of <r>, over the nine MISR viewing angles and 
A ,r pi , A ,r p 2 , 7Vp3, A r p 4 - A> 5 , A' r p6, A >7 and A> x are the numbers of MISR pixels for each 
interval of <r>. At 0° < 0q < 45°, the retrieved r is slightly dependent on satellite 
viewing geometries, particularly when <r> < 10 (see figure 1 2(a)). The agreement 
between the retrievals developed at small solar zenith angles is related to the phase 
function being less sensitive to ice cloud models in the backscattering directions. For 
45° < 0o < 90° and the relative azimuth angle (pot of the MISR Df camera between 
170° and 180°, retrievals developed over small viewing angles agree more with <t>, 
whereas the retrieved ice cloud optical thicknesses tend to be overestimated and under- 
estimated fore and aft of the local vertical directions (see figure 12(6)). The effect can 
be explained by cloud-side illumination and shadowing by 3D clouds (Davies 1984, 
Loeb and Davies 1997, Loeb and Coakley 1998, Iwabuchi and Hayasaka 2002). 

Figure 13 shows the normalized densities of ice cloud OTDs as functions of 0 
retrieved using a water cloud model (figure 13(a)) and the HC model (figure 13(fi)) 
for 263 657 and 244 611 MISR pixels, respectively. Mie theory (Wiscombe 1977) was 
used to compute the scattering properties of water cloud particles, whose refractive 
indices are 1.3244 + 13.58 x 10~ 7 and 1.2901 + i3.94 x 10^ for X = 0.866 |xm and 
X = 2.13 |xm, respectively. Generally, the RMS values of the normalized OTDs in 
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Figure 12. Average ice cloud optical thicknesses retrieved for each of the nine MISR cameras 
(denoted by the subscripts) using the RC model for 1-day matched MODIS and MISR data 
taken over ocean, 2 July 2009: (a) 0° < 6 0 < 45° and ( b ) 45° < 0 0 < 90° and 170° < <p D[ < 180°. 


figure 13 are much larger than those derived for the two granules examined earlier 
(see tables 1 and 2) because the population of optically thick clouds is larger for the 
global data set. The s resulting from retrievals based on the water cloud model is sig- 
nificantly larger than that from using the ice cloud model. In addition, compared to 
figure 1 3(a), the OTDs derived using the HC model are more symmetrical around a 
zero-bias line. These results confirm that the water cloud model is not appropriate 
for retrieving ice cloud optical properties, a conclusion reached from other types of 
analyses (e.g. Minnis et al. 1993, Doutriaux-Boucher et al. 2000). 

Figures 14 and 15 show the normalized densities of ice cloud OTDs retrieved using 
the HC, RC and IC models as functions of 0 . The OTDs for the RC model are smaller 
and more symmetrical around the zero-bias line than those for the other two models. 
It is evident from the close-up plot in figure 15 that more MISR retrievals using the 
RC scattering properties are associated with OTDs in the range of —0.01 to 0.01 . 
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Figure 13. Normalized densities of ice cloud OTDs using (a) water cloud model and ( b ) HC 
model for 1-day matched MODIS and MISR data taken over ocean, 2 July 2009. 


Figures 16 and 17 show the normalized densities of ice cloud OTDs using the 
HP, RP and IP models as functions of 0. Compared to the results in figure 14, 
fewer valid retrievals were obtained using the plate models. None of the plate models 
yield the degree of OTD symmetry about the zero bias that was found for the RC 
model retrievals. Overall, among the six ice crystal habits used for the results in 
figures 14-17, the RC model yields the least angular dependence of the r retrieval 
on viewing angle. This is consistent with the initial test of ice cloud models using the 
two MODIS granules. 

Table 3 lists the average RMS errors of retrieved r for all the ice cloud models tested 
using the 1-day MODIS-MISR data set. The smallest value of s, 0.365, is found for 
the M3 model. Figure 18(a) shows that the normalized densities of M3 OTD produce 
a pattern similar to that found using the RC model (figure 14(6)). Figure 18(6) and 
(c) displays the OTD frequency differences between the M3 and HC results and the 
M3 and RC results as functions of 0. The OTD distribution from M3 is quite different 
from that for the HC retrievals. Conversely, the M3 and RC distributions are very sim- 
ilar as indicated by the frequent occurrence of small differences in figure 18(c). From 
the comparison of OTDs between the M3 and RC results and their small differences 
in e from table 3, it can be concluded that the M3 is minimally better at retrieving an 
angularly consistent r. The asymmetry factors, g, of the ice cloud models are listed in 
table 4. For D e = 50 pm. it can be found that g generally decreases when roughened 
surface and internal air bubbles are accounted for, which is consistent with the results 
reported by Yang et al. (2008b) and Xie et al. (2009). At X = 0.866 and 2.13 pm, e does 
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Scattering angle (°) 


0.000 0.175 0.350 0.525 >0.665 


Figure 14. Same as figure 13, except using (a) HC, ( b ) RC and (c) IC models. 


not have a linear relationship with g, as is evident from tables 3 and 4. However, ice 
cloud models with relatively small g normally have better performance in the retrieval 
of ice cloud properties. 

Figure 19 replots the normalized densities of the M3 OTDs for two different ranges 
of <r>, 0-5 (figures 19(a)) and 5-20 (figure 19(6)). From figure 19, it is evident that 
the M3 model performed well in retrieving ice cloud properties for both optically 
thin and thick clouds. Based on the comparison between the data in figures 19(a) 
and (6), the OTDs of optically thin cirrus clouds are concentrated closer to 0. The 
greater errors for optically thick cirrus cloud retrievals may, in part, be due to a greater 
contribution of 3D effects (Iwabuchi and Hayasaka 2002). 

From these results, it can be concluded that among the tested ice crystal distribu- 
tions, the M3 model is the best model for retrieving cloud optical thickness from solar 
reflectance measurements. However, its accuracy is only slightly better than that found 
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Figure 15. Same as figure 14, except only for normalized OTDs between —0.01 and 0.01. 


for the simpler roughened hexagonal ice column model, RC, and three other habit- 
mixture models. All five represent significant improvements over the use of smooth 
HCs and HPs, roughened plates and IPs. The goal of this study was to determine the 
optimal model for simple shapes and combinations of those simple shapes. As noted 
earlier, many other habits and combinations of habits are possible, and they may yield 
smaller dependencies of the retrieved r on the viewing and illumination geometry. 
However, much additional study is required to examine the sensitivity of the retrieved 
optical thickness to more complex shapes and combinations. 


6. Summary 

In this study, we examined a broad range of ice cloud models that consist of smooth, 
roughened, homogeneous and inhomogeneous hexagonal ice crystals with various 
aspect ratios. The sensitivity of cirrus bulk scattering properties and solar reflectances 
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Same as figure 14, except for using (a) HP, ( b ) RP and (c) IP models. 


to these ice cloud models was investigated using the IGOM and DISORT models. 
At shortwave wavelengths, the smoothed scattering phase function was shown for 
ice cloud models with surface-roughened or inhomogeneous ice crystals. It was also 
demonstrated that cloud bidirectional reflectance is sensitive to the particular ice cloud 
habit with the smoothest reflectance fields resulting from the use of severely roughened 
ice columns. 

The MODIS cloud reflectances from bands 1 and 7, averaged over 1° x 1° 
latitude-longitude regions, are consistent with those in MISR red and NIR bands. 
An algorithm was developed to retrieve ice cloud optical thickness using matched 
MODIS and MISR data. The MODIS data are used to derive ice cloud particle 
size. The retrieved ice cloud particle sizes are then collocated to the MISR pixel 
for cloud optical thickness retrievals at the nine MISR views. Errors in cloud prop- 
erty retrievals are examined when MODIS retrievals are not provided. For the solar 
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Table 3. Ice cloud models and RMS errors of cloud optical 
thickness for the MODIS and MISR granules at a 1-day interval 
on 2 July 2009. 


Ice cloud 
model 

£ 

Ice cloud 
model 

£ 

HC 

0.403 

Ml 

0.367 

RC 

0.368 

M3 

0.365 

IC 

0.389 

M5 

0.366 

HP 

0.375 

M7 

0.369 

RP 

0.484 

M9 

0.388 

IP 

0.394 

Mil 

0.407 
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Figure 18. (a) Normalized density of ice cloud OTDs using the M3 model, ( b ) OTD M i x - 
OTDhc and (c) OTD M i x - OTD RC as functions of scattering angle for 1-day matched MODIS 
and MISR data taken over ocean, 2 July 2009. 


and satellite viewing angles applied most closely to the operational conditions in 
cloud property retrievals, the relative errors in cloud optical thickness are within 
±10% when the ice cloud model consists of only columns. For the plate ice cloud 
model, substantial errors may arise from the absence of MODIS retrievals of cloud 
particle size. Therefore, particle size-dependent biases can be efficiently removed by 
using the MODIS-retrieved particle size in the MISR retrieval of ice cloud optical 
thickness. 

The ice cloud optical thicknesses are retrieved for the nine MISR viewing geometries 
in the overlap regions of MODIS and MISR granules. The differences between cloud 
optical thickness retrieved at the nine views and their mean are used to characterize 
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Table 4. Ice cloud models and their asymmetry factors for D e — 50 |im. 


Ice cloud model 

g (A = 0.866 |xm) 

g (A = 2.13 |xm) 

HC 

0.7938 

0.8452 
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0.7820 

0.8309 

IC 

0.7862 

0.8389 

HP 

0.9172 

0.9352 

RP 

0.8529 

0.8906 

IP 

0.9104 

0.9286 

Ml 

0.7822 

0.8310 

M3 

0.7868 
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M5 
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Figure 19. Normalized density of ice cloud OTDs using the M3 model: (a) 0 < <r> < 5 and 
(b) 5 < <r> < 20. 


the retrieval accuracy. Ice cloud models containing single and mixed ice crystal habits 
are initially examined using MODIS and MISR data at two 5-min intervals. The ice 
cloud model with RCs has the best overall performance in fitting the reflectances with 
the cloud property retrievals. The ice cloud models are further tested using the day- 
time MODIS and MISR data taken over ocean for the entire day of 2 July 2009. The 
retrieved ice cloud optical thicknesses tend to be overestimated and underestimated 
fore and aft of the local vertical directions. Employing mixtures of inhomogeneous 
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and homogeneous ice crystals with surface roughness is found to further reduce the 
uncertainty in the ice cloud retrievals, but only by a small amount. A new ice cloud 
model containing a mixture of ice crystal habits or a distribution of severely rough- 
ened hexagonal ice columns is recommended for future satellite-based retrievals of ice 
cloud properties for both optically thin and thick clouds. 
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